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ABSTRACT 

Aims. We examine the constraints that can currently be obtained on standard cold dark matter models from the most currently used 
data set: CMB anisotropics, type la supernovae and the SDSS luminous red galaxies. We also examine how these constraints are 
widened when the equation of state parameter w and the curvature parameter Q. k are left as free parameters. Finally, we investigate 
the impact on these constraints of a possible form of evolution in SNIa intrinsic luminosity. 
Methods. We obtained our results from MCMC analysis using the full likelihood of each data set. 

Results. For the ACDM model, our 'vanilla' model, cosmological parameters are tightly constrained and consistent with current 
estimates from various methods. When the dark energy parameter w is free we find that the constraints remain mostly unchanged, 
i.e. changes are smaller than the 1 sigma uncertainties. Similarly, relaxing the assumption of a flat universe leads to nearly identical 
constraints on the dark energy density parameter of the universe CI a , baryon density of the universe Sib, the optical depth r, the index 
of the power spectrum of primordial fluctuations n s , with most one sigma uncertainties better than 5%. More significant changes 
appear on other parameters: while preferred values are almost unchanged, uncertainties for the physical dark matter density Cl c h 2 , 
Hubble constant H and <x 8 are typically twice as large. The constraint on the age of the Universe, which is very accurate for the 
vanilla model, is the most degraded. We found that different methodological approaches on large scale structure estimates lead to 
appreciable differences in preferred values and uncertainty widths. We found that possible evolution in SNIa intrinsic luminosity does 
not alter these constraints by much, except for w, for which the uncertainty is twice as large. At the same time, this possible evolution 
is severely constrained. 
Conclusions. We conclude that systematic uncertainties for some estimated quantities are similar or larger than statistical ones. 
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k" ■ 1. Introduction degeneracies. 

Among the many different cosmological tests available 

The field of cosmology has made impressive developments in the beyond the CMB two ^ particu i ar i y robust and based on 

g , past decade. It is now widely accepted that the present universe ge0 metrical considerations: the already mentioned Hubble 

• is well represented by a dark energy dominated universe with diagram of type Ia sup ernovae, and the recent large galaxy 

dark matter being in the form of a cold component and matter surveys that provide d a first detection of the p redicted baryonic 

fluctuations arising from primordial adiabatic Gaussian fluctua- ac0 usticosciilations (lEisenstein et a O005l E05 here after: 

tions generated in the early universe, a picture consistent with |Tegmark et al.^006l T06 hereafter; IPercival et al. 20071 All 

simple inflationary models. these observations reinforce the case of the concordance model 

Such a paradigm was obtained after a remarkable advance in i n cosmology, in which the Universe appears to be in an accel- 

the quality of observational data. The Hubble diagram derived e rated expansion due to some "dark energy" whose true nature 

from type Ia supernovae allowed to obtain the first evidence for remains one of cosmology's most compelling mysteries. The 

dark energy and provided measurements on the proper ties of this ability of the concordance model to accommodate recent high 

component (IPerlmutter et al. 1999I IRiess et al. 1998] ). Another precision data is indeed impressive. Although alternative views 

important contribution came from the resul ts of the Wilkinson exist, it is fair to say that the so called concordance cosmology 

Microwave An isotropy Probe (WMAP) ( JSpergel et al. 2003| has led to specific predictions that were verified a posteriori. 



Spergel et al. 2007) , which have provided strong constraints The most common interpretations for this "dark energy" are an 



on many cosmological parameters and allowed for precise unknown energy component of the Universe or a modification 

estimations of cosmological parameters. Those constraints are i n the equations of general relativity (for some reviews on 

consistent with a wide set of constraints coming from different this subject, see Frieman et al. 2008, Padmanabhan 2006). In 

cosmological tests and have considerably reinforced early claims or der to distinguish between models beyond the ACDM model, 

of flat cosmology (Lineweaver et al. 1993 [Webster etap998). one needs to have very good constraints on the equation of 

The use of combined constraints is therefore a way to increase state of dark energy by combining reliable results from quality 
the accuracy of cosmological parameter estimations and reduce 



L.D. Ferramacho et al.: Cosmological constraints after WMAP5 



observations. Some multi-observational constraints have been 
published recently (Tegmarket al. 2006, Komatsu et al. 2008 , 
Kowalski et al. 2008; [Wright 2007[ |XiaJLaL2QQ8j, but do not 
include a detailed likelihood analysis of the three most recent 
observational data sets (i.e. information from one or several data 
sets was reduced to a single number). 

In this paper, we propose to constrain several cos- 
mological parameters, including the w parameter in the 
equation of state of dark energy (supposed constant), us- 
ing the recent results on the matter power spectrum of 
LRG galaxies, SNIa and WMAP-5 year data. Although 
there are other promising methods to test cosmological 
models and dark energy properties, such as galaxy clus- 
ters dVikhlinin et al. 20"08l lAUen et al. 20081 ) or weak lensing 
(Benjamin et al. 2008 1, we have not included these in our analy- 
sis due to the level of incertitude that may still be present in these 
methods (Ferramacho & Blanchard 2007; Vauclair et al. 2005, 



INuza & Bl anchard 2006; |Blanchard & Douspis 2005] ). Given 



the high number of parameters to be constrained, we adopt the 
Markov chain Monte Carlo (MCMC) method to perform our 
analysis. The paper is structured as follows: In section 2 we 
present and discuss the quality of the data. In Section 3 we de- 
scribe the method and the cosmological models to be constrained 
and in Section 4 we present the obtained results. Our conclusions 
are given in Section 5. 



2. Data 

2.1. CMB 

We use the 5-year data from the WMAP team available at web- 
site http://lambda.gsfc.nasa.gov which consists of the all-sky 
anisotropy maps of the cosmic microwave background (CMB) 
in fits format together with a likelihood code which allows us to 
read these maps and calculate the TT and TE power spectra and 
use them to obtain the likelihood of a given cosmological model. 
We also include the recent ACBAR data dReichardt et al. 20081 l 
which provide additional stringent data on small scales. The new 
generation of codes publicly available that compute very effi- 
ciently the predicted CMB anisotropics in the line of CMB FAST 
(Seljak & Zaldarriaga 1996) have revolutionised the approach of 
constraining cosmological parameters from CMB data. 



2.2. SDSS LRG galaxies 

The shape of the present day matter power spectrum depends 
on conditions that prevailed during inflation (or more gener- 
ally during very early stages of the universe) producing the ini- 
tial fluctuations and on the detailed contents of fields and mat- 
ter of the universe which alter the growth of fluctuations dur- 
ing the expansion. This complex process is summarised through 
the transfer function. In simple inflationary cold dark matter 
models, the initial spectrum is assumed to follow a power law 
: P(k) = Ak" and the final power spectrum can be deter- 
mined using the transfer function which is specified as soon as 
the cosmological parameters are set, including neutrino prop- 
erties. These transfer functions can be obtained through public 



codes like CM BFAST (Seljak & Zald arriaga 1996) or CAMB 
dLewis et a l. 2000). The distribution function of galaxies on 
large scales in the present day universe is supposed to be rep- 
resentative of the present distribution of dark matter providing 
that one can specify a biasing scheme for galaxy formation. The 
SDSS survey dYork et al. 20001 ) allowed to obtain unprecedented 



data on the large scale distribution of galaxies. In this paper 
we use the power spectrum published in T06 and derived from 
58,380 Luminous Red Galaxies (LRG) along with narrow win- 
dow functions and uncorrected minimum-variance errors. For 
comparison purposes, we will also use the data from E05 on the 
correlation function derived from 46,748 LRG, with the covari- 
ance matrix between data points. Both these data sets can thus be 
used for likelihood computation, after providing a model which 
takes into account various corrective effects, such as non-linear 
evolution and redshift-space distortions. 

2.3. SNIa 

Since the first results that showed evidence for an acceler- 
ated expansion of the Universe using a relatively small num- 
ber of distant supernovae, the number of supernovae avail- 
able for determination of cosmological constraints has signifi- 
cantly increased due to several supernova observation programs. 
|Kowalski et al. (2008) have compiled the most up to date set 
of supernovae including the recent large samples of SNIa from 
SNLS dAstieretal.2 006) and ESSENCE (Mikn aitis etal. 20071 ) 
surveys, older data sets and the recently extended data set of dis- 
tant supernovae observed with HST The compilation contains 
414 SNIa and after selection cuts reduces to 307 SNIa. In this 
paper we use the supernovae magnitude and covariance data pro- 
vided by Kowalski et al. (2008) on the 307 SNIa which passed 
selection cuts. 



3. MCMC analysis and results 

The constraints were obtained using the public code 
COSMOMC dLewis & Bridle 20021 ), which implements a 
Metropolis-Hastings algorithm to perform a Markov chain 
Monte Carlo (MCMC) using the above data. This method is the 
one most currently used; we do not address the question of the 
difference between Bayesian and frequentist approaches, which 
could lead to differences in preferred values and uncertainty 
estimates (see for instance Dunkle y~et al. 2008) ). We investigated 
four different sets of parameters to be constrain by this method, 
the simplest being our "vanilla" model with p=(f2 c /i 2 , Q/,/i 2 , 
9, t, n s , A s ). The 6 parameter represents the ratio between the 
sound horizon and the angular diameter distance to the last 
scattering surface. It is used instead of the Hubble parameter 
h since 6 provides a better efficiency and less correlation with 
other parameters. The physical meaning of all constrained 
parameters is described in Table 1 . By imposing a flat Universe 
with w - -1, it is possible to describe all the used data with 
just these six parameters. These parameters were directly 
constrained by COSMOMC, while other important parameters 
as h or Qa were derived from those six basic parameters. The 
two parameters required to correct for non-linearity and bias 
effects (Tegmark et al. 2006} were introduced in the code that 
computes the SDSS likelihood, but marginalised analytically. 
The other parameter sets consisted of adding the dark energy 
equation of state parameter w and the curvature parameter Q.k to 
our vanilla model. 

In order to obtain a better and faster convergence, we ran a 
first chain up to around 10 000 points and then used the results 
to obtain a good estimate of the covariance matrix for the param- 
eters. This matrix was then used as an input to perform a final 
chain run until convergence. 

We ran one chain for each parameter set until it satisfied the 
Raftery-Lewis convergence test (Raftery & Lewis 1992 1 which 
typically happened after 15 000-30 000 accepted iterations. The 
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Fig. 1. Posterior distributions for the parameters constrained in our MCMC analysis. Doted curves represent the mean likelihood of 
the samples, while solid lines show the fully marginalized posterior. The four parameter sets constrained are represented by different 
colours: Vanilla (black), Vanilla + O^(green), Vanilla + w (blue), Vanilla + Q.^ + w (red). 
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Fig. 2. 2D constraints on most correlated cosmological parameters when Vanilla+Q^ + w is considered. The curves show the 
marginalised \cr and 2<x confidence regions 
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Parameter 


Vanilla 


Vanilla + £l k 


Vanilla + w 


Vanilla + fl k + w 


Parameter Definition 


Q. b h l 


0.0227 ± 0.0005 


0.0227 ± 0.0006 


0.0228 ± 0.0006 


0.0227 ± 0.0005 


Normalised baryon matter density times h 1 


Q c h 2 


0.112 ±0.003 


0.109 ±0.005 


0.109 ±0.005 


0.109 ±0.005 


Normalised CDM density times h 2 





1.042 ± 0.003 


1.042 ±0.003 


1.042 ±0.003 


1.042 ±0.003 


Ratio of sound horizon to angular diameter distance 


T 


0.085 ± 0.017 


0.088 ±0.017 


0.087 ± 0.017 


0.088 ± 0.017 


Reionization optical depth 


n s 


0.963 ± 0.012 


0.964 ± 0.013 


0.967 ± 0.014 


0.964 ± 0.014 


Primordial spectral index at k = 0.05Mpc~' 


hg(10 l0 A s ) 


3.07 ± 0.04 


3.06 ± 0.04 


3.06 ± 0.04 


3.06 ± 0.04 


A s is the primordial scalar power at k = 0.05Mpc~' 


n k 





-0.005 ± 0.007 





-0.005 ±0.0121 


Spatial curvature 


w 


-1 


-1 


-0.965 ± 0.056 


-1.003 ±0.102 


Dark energy EoS (w = Pa/Pa) 


n A 


0.738 ±0.015 


0.735 ±0.016 


0.739 ±0.014 


0.733 ± 0.020 


Dark energy density (h z pa) 


Age 


13.7 ±0.1 


13.9 ±0.4 


13.7 ±0.1 


13.9 ±0.6 


Age of the universe in Gyr 


CK 


0.262 ± 0.015 


0.270 ±0.019 


0.261 ±0.020 


0.272 ± 0.029 


Normalised matter density 


C"8 


0.806 ± 0.023 


0.791 ±0.030 


0.816 ±0.014 


0.788 ± 0.042 


Matter density fluctuation amplitude 


<-re 


10.9 ±1.4 


11.0 ±1.5 


11.0 ±1.5 


11.0 ±1.4 


Reionization redshift 


h 


0.716 ±0.014 


0.699 ± 0.028 


0.713 ±0.015 


0.698 ± 0.037 


Hubble parameter 



Table 1. Summary of the posterior distribution mean values for the 
68% confidence intervals. 



different sets of parameters constrained, with the corresponding 



posterior distribution was then used to derive the confidence 
intervals for the parameters. A flat prior was attributed to all 
directly constrained parameters, and there is a restriction on 
h (0.4 < h < 1.0) imposed in CAMB. Table 1 shows the 
marginalised mean values for all parameters constrained, either 
directly or indirectly, as well as the lcr confidence intervals. We 
note the fact that the quoted values do not represent exactly the 
best fit model in either case, but the mean and standard deviation 
of the marginalised probability distribution of each parameter. 
These distributions are shown in Fig.l for comparison between 
the different parameter sets. We also present the 2D marginalised 
distributions for the main cosmological parameters in the case 
where 8 parameters were directly constrained (Fig. [2j». This fig- 
ure allows us to identify possible correlations among parameters. 
Finally, we compared the data published on the evo- 
lution of the Hubble parameter H{z) jSimon et al. 2005 , 



Gaztanaga, Cabre & Hui 2008 1 with the region allowed by our 
constraints on the vanilla model (figure[3]). This was done by tak- 
ing the envelop of H(z) of models falling inside 1 <x confidence 
intervals in the distributions for Q m Qa, Q c /i 2 . This comparison 
illustrates the ability of the ACDM model to fit accurately vari- 
ous type of astrophysical data, although the evaluation of possi- 
ble systematics is difficult for such type of data. 
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Fig. 3. Comparison between the data on the Hubble 
parameter evolution from three data set (triangle: 
Freedman et al.(20 0T)} dark circles |Simonrt"a l. 2005; open 
circles : Gaztanaga, Cab re & Hui 2008j > and one a envelop from 
our constraints on ACDM vanilla model. 



4. Systematics 

4.1. LSS analysis and compressed data 

Table 1 provides an up to date estimation of cosmological 
parameters by combining the likelihood from each data set. 
Clearly, the accuracy is impressive. Given this precision, it is 
vital to evaluate possible systematics that could seriously alter 
these estimations. 

One possible source of systematic uncertainty lies in the 
analysis of the large scale galaxy distribution. Most of recent 
analyses include SDSS data on BAO. In E05 and T06, the au- 
thors have used different descriptors to quantify the galaxy clus- 
tering. In principle, the use of the correlation function or power 
spectrum should give the same results, but using different nu- 
merical estimators can lead to different answers, providing an 
indication on the amplitude of possible biases. Another sensible 
aspect is that of non-linear corrections. In order to compare the 
observed galaxy clustering with that from underlying dark mat- 
ter, one has to account for changes in the linear power spectrum 
(or correlation function) due to non-linear gravitational collapse 
and scale dependent bias of galaxies. Such effects introduce dis- 
tortions on the linear power spectrum shape and can even shift 
the position of the acoustic peak of up to 3% (Smith et al. 2008). 
Several models have been proposed to correct for these effects, 
usually supported by numerical simulations. In T06, non-linear 
corrections were performed using a one parameter model intro- 
duced bv lColeetal. 20051 



* gal — *dewigg& 



1 + Q„,k 2 
1 + 1.4k 



(1) 



Here Pdewigg stands for the linear power spectrum corrected 
for a power suppression of the baryonic oscillations due to the 
displacement of galaxies induced by their particular velocities. 
The other part of the above expression quantifies the scale de- 
pendent distortion due to non-linear gravity and bias. 

In E05, the authors used the same correction for the bary- 
onic peak supression, but the scale dependent bias and non-linear 
collapse were treated differently, by means of the HALO FIT for- 
malism and an expression derived from the numerical simulation 
of |Seo et al. (20 08) and applied directly to the real space corre- 
lation function. 

We have thus compared the different approaches referred to 
above, when combined with the CMB and SNIa information. 
We restricted this analysis to the vanilla model and the results 
are shown in Fig. 3 for relevant parameters. One can see that 
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Fig. 4. Joint marginalised posterior distributions for relevant cosmological parameters in a vanilla model comparing four approaches 
to the LSS data: Matter power spectrum and corrections for non-linear effects from T06 (black), correlation function from E05 with 
corrections for non-linear effects from E05 (red) and T06 (blue). The final constraint uses the distance parameter A derived in E05 
(green). 



using the E05 two point correlation function and non-linear cor- 
rections results in a difference in the mean value of some param- 
eters compared to using the T06 data on the matter power spec- 
trum together with a Q-model. The parameters most affected are 
the densities Q c h 2 , O m , £2a and the Hubble parameter Hq. For in- 
stance, Q.„, is shifted from 0.262 to 0.275 while the uncertainty 
is of the order of 0.015, the difference is therefore of the order of 
1 <t or more. 

In order to track the origin of such a difference, we ran an- 
other chain using the correlation function data from E05, but 
with a Q-model correction for non-linear effects. One could ex- 
pect that by using the same modelling as T06, the observed 
difference would be reduced. However, the blue curves in Fig. 
3 show that the difference increases in this case, with 0,„ = 
0.295 + .015. This result is the same if we leave the Q„i pa- 
rameter in eqQ] as a free parameter or if we use the best value 
found by T06, Q„i = 31. Such a result indicates a significant 
difference between the estimations in E05 and T06, since they 
both used similar galaxy catalogs. Although the observed differ- 
ences are no more than 1 to 2 <x, they are significative enough to 
consider that the present LSS and BAO data has to be used with 
some caution. Clearly, future works on this subject must address 
carefully the issue of bias in estimators (Norberg et al. 2008 ) and 
modelling of non-linear effects. With this caveat, we stress that 
the use of the matter power spectrum was less time consuming 
and we kept this approach for wider parameter investigations. 

It is also interesting to notice that several recent papers pro- 
vide constraints from combination of these three data sets, but in 
which some of the constraints were obtained from highly com- 
pressed information, typically one number, for instance in the 
form of a reduced distance. Obviously different compressed in- 



formation from the same data could result in different constraints 
(see for instance Lewis, 2008); the use of compressed informa- 
tion and approximate likelihoods should therefore be tested to 
ensure they do not lead to inadequate approximation of the like- 
lihood on the full data set. Although some compression is nec- 
essary to reduce the numerical cost of likelihood analyses, too 
strongly compressed information may result in biases on both 
the preferred values and on the uncertainties. We have used the 
distance parameter A from E05 to perform our constraints, with 
the results also plotted in Fig. 3 (red curves). This parameter 
seems to capture quite accurately the cosmological information 
from the full correlation function with non-linear corrections 
used in E05. The mean values for all parameters are highly com- 
patible in both analyses. The main difference lies in the uncer- 
tainties of some parameters, which seem to be overestimated by 
an factor of 25% to 30% when we use this reduced distance. 



4.2. Possible astrophysical troubles 

Even in the presence of ideal data with no systematics or in the 
data acquisition either in the data analysis, inferred constraints 
may be biased because of unidentified complex astrophysics. For 
instance CMB fluctuations are contaminated by various fore- 
grounds, which may introduce bias in estimation if incorrectly 



subtracted (Douspis et al. 20061. Similarly, the acoustic scale as 
well as the whole correlation function inferred from galaxy sam- 
ples could be biased by distortion in redshift space, non-linear 
dynamics and the biasing of the galaxy population. All these ef- 
fects are believed to be controlled to a much higher precision 
than present surveys allow by now (Seo et al, 2008). 
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4.2.1. SNIa evolution and consequences 

The situation for SNIa is more problematic: SNIa are com- 
plex objects probably resulting from mass accretion of a white 
dwarf in a binary system; the properties of the final explosion 
may depend on many of the parameters of the progenitors (age, 
metalicities, environment, ...) and the properties of the explo- 
sion are not well known (detonation/deflagration, physics of 
the ignition points). The SNIa results heavily rely on the as- 
sumption that the peak luminosity, once corrected for duration- 
luminosity dependence (Phillips 1993), does not vary with red- 
shift. Some evolution in the progenitor population is expected 
from the general stellar population. Some evolution may there- 
fore be present and hard to detect. Indeed, this is one of the main 
weakness of SNIa results dRiess & Liv io 2006). From the ob- 
servational point of view only differences in the spectroscopic 
properties of distant and close supernovae can be quantified 
(Balland et al, 2007; |Foley et al. 2007) . The role of possible 
grey dust (Goobar et al. 2002) or evolving dust may also be im- 
portant (Holwerda, 2008). 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 



Fig. 5. Constraints in the Qa - b (top) and Q.^ - K (bottom) 
planes from the Hubble diagram for a flat cosmological model. 
Contours correspond to 1 and 2 cr regions on two parameters. 



Despite the remarkable efforts on observations of distant 
SNIa, it is clear that the hypothesis of negligible evolution is ex- 
tremely difficult to prove. As it is commonly written, the absence 
of evidence is not an evidence of absence. A simple evolution 
law could be: 



dNordin et al. 2008). Wright (2002) showed that an evolution 
model in which peak luminosity varies as an exponential func- 
tion of cosmic time within an Einstein de Sitter universe may 
mimic the accelerated expansion. We have therefore also exam- 
ined the possible consequences on cosmological constraints as- 
suming that the peak magnitude of SNIa evolved linearly with 
time as: 



Am(z) = K 



to - Kz) 
t -h 



(3) 



Am(z) = bz 



(2) 



where to is the present age of the universe, t\ is the age of the 
universe at z — 1 and t(z) is the age of the universe at the redshift 
of the supernova explosion. The parameters b or K represent the 
typical change in magnitude at redshift 1. We first examined the 
consequence of such evolution for the SNIa Hubble diagram. 
In Figure [5] we determined the constraints in the O a - b and 
Qa - K planes from the Hubble diagram for a flat cosmolog- 
ical model. As one can see, there is a degeneracy between the 
two parameters and no tight constraint can be obtained on Q. A 
any longer. As the second evolution law introduced stronger lev- 
els of degeneracies, we concentrated on this second form. We 
then re-examined the various constraints established in section 
3 with our additional parameter describing SNIa evolution. This 
was first done on the vanilla model with K as an additional pa- 
rameter. We found that the cosmological constraints are almost 
unchanged. This is not surprising as we checked that removing 
SNIa does not change the constraint on the vanilla model by 
much. Our analysis allows us to constraint the amount of possi- 
ble evolution in SNIa in our model consistent with existing data 
(CMB and correlation function). We found that K is essentially 
uncorrected with the other parameters. The posterior distribu- 
tion on K is given in figure [6] for our parameter sets. The con- 
straint on K can be summarised by: 

K = -0.042 + 0.042 

for the vanilla model. We also examine how constraints are re- 
laxed when the curvature parameter, Oa and w are free. Again, 
the constraints are not modified significantly. Also, there is a 
clear degeneracy between w and K, which degrades constraints 
on both parameters. Resulting constraints are summarised in 
Table 2. For the most general model (curvature and w free) the 
constraint on K is three times larger: 

# = -0.133 + 0.077 

Our analysis shows that evolution is not required by the data. 
Most noticeably, it changes significantly the constraints on w; 
while the vanilla model is almost at the center of the 1 <x contour 
when evolution is not allowed (Fig. |8j, the constraints appear 
twice as bad (darker contours in figure [8} when evolution is al- 
lowed, and the best model then lies close to the boundary of the 
2 cr contour without evolution. 

Interestingly enough, constraints on others parameters are 
quite stable even in the most general case (curvature and w free) 
and therefore insensitive to such a possible evolution. Our re- 
sults put a stringent constraint on the possible evolution of SNIa 
and thereby provide an idea of the level at which systematic un- 
certainty for possible evolution should be controlled in order to 
provide 

useful constraints on the problem of the determination of cos- 
mological parameters in conjunction with other methods: even 
in the worse case the amount of evolution is found to be of the 

order of Am 0.136(1 cr level), at redshift z - 0.5, while in the 

vanilla model we found that evolution is limited to Am 0.055 
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Parameter 


Vanilla 


Vanilla + il k 


Vanilla + w 


Vanilla + £2 k + w 


Qth 1 


0.0228 ± 0.0006 


0.0227 ± 0.0005 


0.0227 ± 0.0006 


0.0226 ± 0.0006 


Q, c h 2 


0.1 10 ±0.004 


0.109 ± 0.005 


0.113 ±0.005 


0.111 ±0.005 


e 


1.042 ±0.003 


1.042 ±0.003 


1.042 ±0.003 


1.042 ±0.003 


T 


0.088 ±0.017 


0.087 ± 0.017 


0.085 ± 0.017 


0.085 ± 0.016 


n s 


0.968 ±0.013 


0.965 ±0.013 


0.963 ± 0.014 


0.960 ± 0.014 


log(W l0 A s ) 


3.07 ± 0.04 


3.06 ± 0.04 


3.07 ± 0.04 


3.06 ± 0.04 


fit 





-0.002 ± 0.007 





-0.017 ± 0.013 


w 


-1 


-1 


-1.112 ±0.148 


-1.33 ±0.242 


K 


-0.042 ± 0.042 


-0.035 ± 0.042 


-0.105 ±0.091 


-0.133 ±0.077 


n A 


0.747 ± 0.017 


0.745 ± 0.020 


0.756 ± 0.022 


0.744 ± 0.022 


Age 


13.6 ±0.1 


13.7 ±0.4 


13.6 ±0.1 


14.5 ± 0.7 


"m 


0.253 ± 0.017 


0.257 ± 0.025 


0.244 ± 0.022 


0.272 ± 0.029 


G"8 


0.801 ± 0.026 


0.794 ± 0.029 


0.846 ± 0.068 


0.867 ± 0.060 


<-re 


11.1 ±1.5 


11.0 ±1.4 


10.9 ±1.5 


10.8 ±1.4 


h 


0.725 ±0.017 


0.720 ± 0.036 


0.748 ± 0.038 


0.703 ± 0.042 



Table 2. Similar to Table 1, but including an evolution parameter for the supernova intrinsic luminosity. 




-0.5 



Fig. 6. Posterior distribution for the SN evolution parameter K, 
for different parameter sets of CDM models. The colour code is 
the same as the one used in Fig. 1 . Evolution is systematically 
preferred, although the non evolving solution remains accept- 
able. 
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0.02 0.04 



Fig. 7. Joint 2D marginalised constraint on the dark energy equa- 
tion of state parameter w (supposed constant) and the curvature 
density Q.^, assuming evolution in the SNIa luminosity (blue) 
and considering no evolution (red). The contours show the con- 
fidence regions of 1 (darker) and 2 <x. 



at the same redshift. We notice that spectral SNIa evolution is 
constrained to be less than 10% from the observed properties of a 
sample of SNIa with a median redshift of 0.5 (Foley et al. 2007). 
Finally, we have examined how the so-called figure of merit 
(f.o.m), defined according to the DETF (2006) (see for instance 
Virey & Ealet 2007), is modified when possible evolution in the 
SNIa of the above type is taken into account. We therefore used 
the CPL formulation (IChevallier & Polarski 20"0T1 Hinder 20031 
for dark energy evolution. The state equation parameter is sup- 
posed to evolve according to: 



w(z) = wo + (a - \)w\ 



(4) 



where a is the expansion factor normalised to one at the present 
epoch. We restricted ourselves to the fiat case. The f.o.m is de- 
fined as the area of the 90% contour in the wo, w\ diagram. The 
joint constraint contours are presented in figure [8] The f.o.m. is 
found to be 2.03 when no evolution was assumed and is low- 
ered to 1.17 when evolution is considered. This is in agreement 
with our previous finding : taking into account evolution leads to 
degraded constraints although the loss in precision is limited. 



S 
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Fig. 8. Joint 2D marginalised constraint on the dark energy equa- 
tion of state parameter w (supposed constant) 
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5. Discussion and conclusions 

In this paper, we have used three of the most solid sets of obser- 
vational data in cosmology to perform an up to date constraint 
on the ACDM model: CMB, the Hubble diagram of distant SN 
and the shape of the power spectrum of galaxy distribution on 
large scales. We have chosen to use the power spectrum shape 
because the analysis using the correlation function was apprecia- 
bly longer. We have confirmed that 6 parameters only are enough 
to reproduce the three currently used data sets. Using appropri- 
ate likelihoods for each data set, we obtained very good con- 
straints with typical uncertainties of the order of a few percent. 
The comparison between the values obtained in this work and 
those found using only WMAP3+SDSS ( Teg mark et al. 2006) 
data do not show significant changes in most of the parame- 
ters, which demonstrates the great constraining power of these 
large scale observations. Indeed, the introduction of the super- 
nova data produces changes mostly on the parameters that gov- 
ern the geometry of the Universe, namely Q.^ and h, which have 
their mean values reduced by around 0.02/0.03. As for the dark 
energy EoS parameter w, we do not obtain a better constraint by 
adding the supernova data, and the results point to a value highly 
compatible with a cosmological constant although the \cr error 
bars are appreciably larger than on other parameters (of the order 
of5tol0%). 

We found that adding extra freedom to the vanilla model 
does not decrease the quality of the constraints in an apprecia- 
ble way for several parameters: as can be seen from figure Q] 
constraints on the baryonic content of the universe Q/,, the dark 
energy content Qa, the matter content of the universe Q. c and the 
primordial index of scalar fluctuations ns are almost unchanged. 
We get a slightly wider distribution in other parameters, with 
preferred values unchanged but with uncertainties increased by 
a factor of up to two. The most drastic change happens for the 
age of the universe, with an uncertainty increased by a factor of 
more than five, although the uncertainty remains still only of the 
order of 5%. 

Our constraints have been established within some assump- 
tions. For instance we do not investigate the possible contribu- 
tion of tensor modes (Zaldarriaga & Seljak 1997; Liddle & Lyth 
2000) or neutrino mass (Goo bar et al. 20 06) . Obviously, the ad- 
dition of new freedom could modify some of our constraints. 
There is no general method to avoid this type of limitation; for 
instance, one of the most critical assumptions on constraints ob- 
tained from the CMB is the power law shape of the spectrum of 
primordial fluctuations. Relaxing this assumption could lead to 
dramatic modification of the constraints (Blanc hard et al. 2003, 
IHunt & Sarkar 2007) . 

While preparing this paper, we benefited from the release of 
the 5-year WMAP results. The W MAP team presents M CMC 
constrains on the ACDM model (Komatsu et al. 2008) . using 
WMAP+SN+BAO They used the BAO measurements on the 
angular diameter scale Da derived from the SDSS LRG corre- 
lation function (Eisenstein et al. 2005), while we preferentially 
used the full power spectrum shape data. By comparing the two 
approaches we find very similar constraints on all parameters 
and similar l<x confidence intervals, which indicates that the use 
of the distance parameter A in implementing the LSS constraint 
is essentially sufficient. Also, we find that the use of the power 
spectrum instead of the correlation function leads to appreciable 
difference of up to two sigma. Finally our investigation on the 
consequence of a specific model for the evolution of SNIa lumi- 
nosities has been interesting in several respects: first, although 
evolution does not allow any conclusion on cosmological pa- 



rameters from the Hubble diagram alone, we found that most of 
the constraints from combinations remained unchanged, and at 
the same time the amount of evolution is severely constrained to 
be less than 5 to 20% for redshifts between 0.5 to 1 . 

The current status of cosmology seems then to be highly 
compatible with the ACDM. Present-day existing data allows 
to estimate cosmological parameters to within a few percent 
whithin a specific model. This raises the fundamental question 
of the role of astrophysical sources of uncertainties. Clearly the 
next step is to understand in more detail the nature of dark en- 
ergy (cosmological constant, dynamic scalar field or modified 
gravity) and this is yet to be achieved. The ever increasing qual- 
ity of data, anticipated from future large observational programs, 
could be used to rule out or validate some of theses alternatives. 
However, future projects would need to keep systematic effects, 
both from methodological and astrophysical origin, controlled to 
achieve a typical anticipated precision of one percent or better. 
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